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© A system for predicting the life of a cutting toolln.afljiujgrrxa ted metal cutting machine, com prises a plurality 
of sensors for measunr rgzM ysical data (601) re lating to a cutting operation of said machine, and a comptrter-- 
connected to "these sensors. The computer comprises a feature extractor (603) for transforming the physical data 
into feature values, and a tool life predictor (605) for arriving at a prediction of the Jife^qf the cutting tool. The tool 
life predictor has an influence diagram having an input level of feature nodes^Jntermed^ 
the feature nodes and to a machining class node t -and-an-ootplKnn^e^6?mected to the intermediate nodes. The^" 
tool life predicator further comprises a first module for computing averages of feature values for each sensor, a 
second module responsive to the feature averages for making a classification of the feature average as indicative 
of a long life or a short life, and a third module connected to the influence diagram for combining the 
classifications to arrive at a classification and a conditional probability at each intermediate node and a 
classification at the output node which is the predicted life of the cutting tool. 

Furthermore, a method to enable a real-time expert system to learn is provided which comprises the steps of 
creating a relational structure corresponding to the influences that a first and a second variable have on a third 
variable, obtaining data for these variables over some duration- of time, determining a tolerance zone around a 
sample which divides classes and using a dynamic distance classifier, which, is responsive-to the tolerance zone, 
to determine decision boundariesTor the first and* second variables "and. probabilities associated with the decision 
boundaries. The probabilities indicate the confidence with which a .decision boundary classifies values of- the first 
and second variables with respect ttrthe va|ue of the* third variable. * ■ - - 
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NOTICE * " 

(C) Copyright- Texas Instruments Incorporated 1991. A portion of the disclosure of this patent document 
contains material which is subject to copyright protection. The copyright owner has no objection to the 
5 facsimile reproduction by anyone of the patent document or the patent disclosure, as it appears in the 
Patent and Trademark Office patent file or records, but otherwise reserves all-eopyrights whatsoever. 

BACKGROUND OF THE INVENTION ~ * T ' 

70 1 . Held of the Invention 

This invention relates to manufacturing automation and in particular to a system and method for timely 
replacement of cutting tools through early prediction of cutting tool life and accurate diagnosis of tool wear 
using an Influence Diagram Expert System and multiple sensors. 
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2. Description of the Related Art 



The background of this invention is described in conjunction with manufacturing automation as an 
example. iv 

20 Automation in manufacturing is a rapidly growing field with significant impact on improved product 

quality, high productivity and reduced manufacturing cost. While it has been adopted successfully at 
different levels of manufacturing process, one of the main hindrances of completely^utomating any metal 
cutting operation is the timely replacement of cutting; tools. Currently, hurrVan; intervention is required to 
change the cutting tools at appropriate heuristicany chosen fntervals and "typically, one machinist is 

25 assigned to supervise and replace cutting tools fbr'a batch of automated machining centersTWith adequate 
sensors for in-process monitoring of the tool condition, ah, improvement can be obtained in machining 
economics, especially in high volume production lines. Cutting tool condition monitoring r wear diagnostics 
and appropriate machine control problems have been investigated through different methodologies: real 
time expert system based approaches and mathematical model based approaches. While these on-line 

30 monitoring techniques provide a solution to correctly identify a damaged or worn tool, they are postmortem 
techniques; waiting for an appropriate decision after the tool is worn or damaged. 

A tool changing operation may involve multiple tasks like procurement of a replacement tool from a 
centralized storage area, loading tool-setting software from the plant manager to the machine controller, 
dynamic re-scheduling of the production line to account for the anticipated machine down time, and so on. 

35 The present tool monitoring techniques do not provide adequate lead-time to initiate these multiple tasks in 
preparation for a tool-change, and hence, no significant improvement is made in reducing the downtime due 
to cutting tool changes. Another added limitation of present continuous on-line monitoring is the requirement 
of dedicated processors on each machine, primarily because of stringent response time requirements of the 
present day real time diagnostic systems. 

40 In order to automate the cutting tool replacement task with adequate forewarning, a reasonably accurate 

estimate of the life of the cutting tool is required. As mentioned earlier, a significant advantage would be 
that the predicted tool life information could be utilized for dynamically re-scheduling machining operations 
and cutting tool replacements with a reduction in inventory and labor costs. 

Mathematical model based techniques, called Taylor's tool life equation, and extensions thereto, may 

45 be applied to estimate the cutting tool life. However, these equations provide a poor estimate of the life of 
the cutting tool. For example, in the case of drilling, work piece material hardness is an important factor 
affecting the life of the drill but is not taken into account. While modified forms of Taylors equation with 
corrections for work piece hardness have been applied to determine the drill tool life, such an approach is 
impractical in a production line since hardness for each work piece is not readily available or easily 

so measurable. Furthermore, in some other machining situations, Taylor's model has been found unsuitable. 

Therefore, it would be desirable to have a system in which a cutting tool's life can be accurately 
predicted during the initial use of the tool so as to provide adequate lead-time for the replacement of the 
tool. 

In the area of cutting tool wear diagnosis, some pr s nt day techniques are pattern recognition (PR), 
55 n ural nets (NN) and real-time xpert systems (RTES). PR and NN, use pattern sampl s that are already 
classified based on some wear criteria and train the system to recognize this fact. This approach is called 
'sup rvised learning'. In 'unsupervised learning*, no wear crit ria is used and hence, the knowledge about 
the class to which the samples belong is not available a priori. Such an unsupervised training approach 
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using fuzzy clustering, technique has been applied to the drill. tool classification problem. Using.; thrust and 
torque data from a torque dynamometer, a fuzzy classification of the state , of, -the. drill tool has been 
performed successfully. .However,, the drill is classified as. w.grn only after the cutting, e^dges are severely- 
damaged. This prediction is delayed until the dri II is about to fail. . v. ~ . 

5 The advantages of the RTES techniques, are the. ease of modification, the,. ability t to -select optimal 

control decisions jby- optimizing over, a cost function after classifying -the current state.. of the tool, and the, 
facility to repetitively utiljze the.- same inference engine : ; while applying, to different machining operations, If. 
the diagnostic system needs to be updated to. handle, new combinations of machining^ parameters pr^new, 
sensors, then in case of PR or NN techniques, the system has : to rbe retrained with new as well as old data, 

jo a rather time consuming process. However, in case of RTES, the knowledge base can be. easily updated 
with the new information. - ■ r .. 

Heretofore, one of the problems with the present RTES technique has been the. lack of learning . 
capability. The; development and tuning of ther.RTES requires detailed analysis of* pertinent data to extract, 
features and subjective f -innate of the conditional-probability distributions that form the critical elements of. 

75 the knowledgebase.; , 

SUMMARY OB THE INVENTION . t , cr: V 

A system for predicting the life of a cutting tool in an automated metal cutting machine, comprises a 

20 plurality of sensors for measuring physical data relating to a cutting operation of said machine, and a 
computer connected to these sensors. The computer comprises a feature extractor for transforming the 
physical, data into feature values, a tool life predictor for arriving at a prediction of the life of the cutting tool. 
The tool life predictor ,has an influence diagram having . an input level of .feature nodes, intermediate nodus 
connected to-: the feature nodes and to a machining class node, ap output node . connected to ti-^. ; 

25 intermediate nodes, a first module for computing averages of feature values for each sensor,. a secqtyl 
module responsive jo the, feature averages for, making a , classification of the feature average as indicative of . 
a long life or a: < short life, and ;a - third module, connected to the influence .diagram for combining the 
classifications to arrive: at a classification and a conditional probability, at each, intermediate node and a 
classification at the' output node, which- istthe.predicted life of ther.cutting,tooL - 

30 Furthermore, a method to enable a realrtime expert, system, to learn comprises the steps of creating a 
relational structure corresponding-, to the; influences that a :first £ and ^ a - second variable have on a third .- 
variable, obtaining, data for these variables p.verr sorrje durationof time, determining a tolerance zone around . 
a sample which 1 divides .classes ,;and.,- using a dynamic, distance; classifier, :which is responsive to, the 
tolerance zone, .to determine decision boundaries J or the first; and -second, variables and .probabilities,, 

35 associated with the decision boundaries. The probabilities indicating the confidence with which^a decision 
boundary classifies values of .the first and second -variables with respect to the value of the third, variable. 

BRIEF DESCRIPTION} .OF THE DRAWINGS 

40 In the drawings: ; = - - 

FIG 1. is a . block, diagram of a computer aided manufacturing (CAM) system; 

FIG 2. is a* schematic, diagram showing the logical links between sensor data in a CAM system. and an 
Influence Diagram- Expect System (IDES).,,' - <- b *: - ■ 

FIG 3. is a flow charHor operating *a CAM, system using an IDES to predict tool life and to monitor tool 
45 wear; c :\. .. • ; ; ;*•;/, - .• ,■• . .. - _ : . - r _ t ... 

FIG 4. shows the- relationship* of the traini.ng phase r to the operations phase;. 
FIG 5. shows an-exampie influence.diagram; . w . - ■ ■ - ■ •: • - 

. FIG 0. shows the influence., diagram, topology for predicting; tool life according to the present in ventipn; . . 
FIG ? is a graph of experimental data used during the training phase of. the present invention; 
so FIG b. shows the influence diagram topology for monitoring tool wear according to the^prese.nt invention; 

FIG 9. is a flow chart .showing how. a -dynamic distance ^ classifier^ according to the present invention is ? 
utilized in a computerized machining operation, i r -\ . <r. . :.6 - *;;v, ;r 
Corresponding- numerals , and r symbols in the different figures refer ..to-- corresponding,, parts, unless, 
otherwise indicated. . ^ ... , t 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS ..^ vprr ^, ^ =r- 

In Figure 1, is shown- iruan. automated ^manuf aeturing operation using this jnvention. A . metal cutting 
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machine 90 hasia drill bit' 100 attached to-a quill* 108, rotated by a spindle motor 106.' The spindle motor 
106 is connected to ah AC electrical power source 117 via : wires-i 16. 

Drill bit 100 isjowered by feed- motor 110 trkd a work' piece 102 attached to a work table ^26 by vise 
104. Feed motor 110 is connected to a DC electrical power source' 1 19 via wires 120. 

5 A number of sensors measure physical data pertinent to the drilling operation of machine 90. A spindle 

motor current gage 114 senses the AC current transmitted through wires 116, a feed motor current gage 
118 measures 'the DC cuttent transmitted through wires 120, a strain gage 122 -is attached to* quill 108 for 
measuring the thrust force exerted by the tool 100 on work piece 102, an accelerometer ~'T28 senses 
vibrations in vise i 04, and a dynamometer 1 24 measures the force exerted by vise 104 on table 126. The 

w outputs from these sensors-are fed into an analog signal conditioner 130 for adjusting the raw signals into a 
range acceptable as input by analog-to-digital converter 132. Conditioned signals from the analog signal 
conditioner 130 are transmitted through lines 160 through 170 into analog-to-digital converter 132 of a; 
digital computer 134. Computer 134 contains a central processing unit 136 and a memory 138. Connected* 
to computer 134 are monitor 144* keyboard 148, disk drive 140, and tape drive 142. 

75 Adjacent to machine 90 is a measurescope 148 which has a monitor 150 attached to it.: The 
measurescope 148 is positioned so that it can visually measure the wear scars of drill bit 100. 

Computer 134 processes inputs from sensors 114, 118, 122, 124, and 128 to predict when drill bit 100 
will be worn out and at the end of predicted life to monitor drill bit 100 for wear on a continuous basis. The 
prediction and diagnosis of wear may either be displayed on monitor 144, for a human operator to take 

20 some action, such as replacing the tool, or may initiate the automatic procurement and replacement of drill 
bit 100. The replacement may be done by some process control hardware "11 2 attached to computer 134. 

During certain phases of the operation of machine 90, computer 134 will be idle^with respect to 
processing data and controlling' machine 90. That idle time may be used to service other machines^ (not 
shown). Thus,' computer 134 may be used to predict tool life and diagnose -fool wear in multiple rnachines .- 

25 which operate'' irVparallel. • :;: Z/ * > % "- r - ; • . 1 - • ' ^ -^-^ 

The operatiorTof machine 90 is controlled by an Influence- Diagram Expert System *"<IDES) 154 ; -shown in: 
FIG 2. As Figure 2 shows, computer 134 executes Influence Diagram Expert System (IDES) 154.. It predicts: 
the length of life of drill bit 100, monitors drill bit 130 for wear, and initiates the : ^ replacement procedure^ 
when drill bit 100 has been diagnosed as wortf-out" Figured shows a schematic" "of the logical links between 

30 sensor data 156'and IDES 154. The sensor data is,- 'for example,-sensors"TT4, 118; 122, 124/and 128, of -~ 
Figure 1. IDES 1 54~is alsb connected to knowledgebase : 1 52-, in which" IDES 154 stores information it uses- 
to process data obtained from the sensors to mane decisions' which" are displayed- on monitdr 1"44>orwhich 
initiate some action on the part of process control hardware 1 T2T Knowledge base 152 may be physically 
located in memory' 138 Fig. 1, on disk '14d,*oh tape 142;* or distributed -over these devices or other 

35 storage devices not' shown. " y J -.Si.*.. -n. ■ ■ ' : 

IDES 154 is also logically connected to parameters 1 56/ These - m ay r for examplev be entered from . 
keyboard 148, of Fig. i, or retrieved from data files on disk 140 or tape 142. Parameters 156 includes, for 
example, information for categorizing the machining operation, such as, tool size, worlc-piece material, drill, 
feed-rate, and drill speed. 

40 The operation of machine 90 and IDES 154 consists of the following phases: ' * 

1. An experimental phase in which the drilling 'iteration is carried out and sensor measurements are - 
taken at regular intervals throughout the life of the tool. Furthermore, the point- at which the drill has been* 
worn out, i.e., its life has been expended, is noted. As each hole is drilled Ks dumber is Voted and each 
data point is indexed by it. Similarly, the point at which the tool is deemed'to have been worn : out: is also 

45 noted by hole number. This experimental data is later used to build knowledge base 152 for both tool life 
prediction and tool wear diagnosis. The knowledge base 152 may be stored orTdisk 140 or tape* 142: . 

Both raw sensor data and quantities computed from the raw sensor ■ data ate used' in both tool life 
prediction and tool wear diagnosis. The'manner fn 'which the computed quantities' are" derived .from the 
raw sensor data is discussed in 'greater "detail below. The raw sensor data and the computed quantities: 

so are collectively called tea tufes. : " ' i6Ty ' r 

2. A training phase, in which'the knowledge base is : built and ■ fine-tuned cohsists'of two d istincf Aspects: 
first, building the knowledge base for accurate prediction of tool life and ; -second, : the"knowledge base for 
the diagnosis of tool wear. Certain portions £f the^trainirig^phase are automated -by using the Dynamic 
Distanc Classifier algorithm d scribed below. 

55 3. An operations phase in which IDES 154. utilizing the knowledge base which was built during the 
training phase, predicts tool life and monitors tool wear. ' ,l '-"' , " Jv * - • *- • :i 

Thus, th goal is to build a knowledge base which an IDES 154 uses to predict tool life and to monitor 
drill wear. The training phase is mor* asily understood if the operations phase is considered first. 
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OPERATIONS PHASE 

FIG 3. is a flow jchart Rowing the operation; of a computer aided; manufacturing .operation according, to 
the present invention.. A^drilling operation starts* at Step 601. In .conjunction with the:drilljng^ measurements 

5 are taken from the various sensors* used in the Influence DiagranrV to~predict'tOoNife: Early in the-operation, - 
for example after three holes, the features are computed, Step 603. Those features aroused by IDES 1541 

to predict the life of -drill 100, Step 605 _i. * : . 

If the prediction is that the tool will have* a hong life j(Step 607), the-machine r continues drilling until the 
tool's predicted ;( mid-life, Step 609. At which point, a verification . procedure is" invoked x Jf : the earlier' 

io prediction of a I6ng life is verified the~drilling "continues until the end" of the predicted life-oMhe tool, Step 
615. At that poirrt sensor- data 'is - coElected and; features computed* fdrHhe next three holes, Step 619. The 

process cycles bask Jo predicting.^ 

During the time between -prediction of a "long life" (Steps .605 and 607) and the verification at mid-life 
(Step 609) computer 134 does not process data from machine 90,, Thus,. during that time, computer 134 is 

75 available for othertasks, forexample, servicing' other machines attached to itr - — 

If a long life!;is not predicted -in Step 605, or if -a "predicted long life is not verified: at the predicted mid- 
life in Step 609Mhen a replacement tool is obtained, Step 61 3J- One -^embodiment is represented by the 
dashed path from Step, 613JO Step 623. Drilling continues until the predicted end of the life'oftne tool, Step 
625. The tool is replaced at the 'point at which~JDES ~154 has predicted 't|iat:the tool's life erVds, Step 625. 

20 An alternative embodiment is~ indicated by * the- solid path emanating from Step 613. Because -it is 

desirable to never use a -tool beyond the point" at which it is worn out,* thejtpol life prediction^ adjusted to 
be a conservative^ejstimate. _For that reasojr, tqols_ which have-not worn* out may have been predicted as 
worn out. Rathe* than -discarding such tools, .the embodiment represented by the path through Steps 6i3 K 
617, and 619. provides a' method for using the "sensors 114, 1TB', 122,126; 'and 128, to : diagnose the tool for 

25 wear. Thus, IDES^54-eonfcinueusly- monitors-the-sensGrs^arid^^ 

continues, Step 621, until IDES 154 diagnoses .the end of the life of the tool, Step 621. When the tool has 
been expended the tool is replaced, Step 623. The entire' bperatioVi starts over with the freSh drill at Step 
601. 

The method described above in conjunction with Fig. ^3 may be implemented to be executed on 
30 computer system 134. One ; embodiment has a functional module to compute features from the raw sensor, 
data (Step 603), a functional module each for the tool life prediction (Step 605), the verification of the ^ 
prediction (Step 609), and the tool wear diagnostics (Steps 61 7 and 621 ). 
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TRAINING PHASE 



Fig. 4 shows the relationship of the training phase , for tool life prediction to .the operation of the IDES' 
154 and the automated manufacturing system shpwn in Fig.,1. A.-similar approach is used' for training the 
system for tool wear diagnosis. - 

The training operation may be-divided functionally jnto two^major categories:, offline knowledge base 

40 generation and on-line tool-life prediction. The off-iine task generates the qualitative, and the quantitative 
knowledge. base for the influence diagram model of the drilling process. The knowledge base is used by the 
real-time system for on-line.-, prediction of tool-life and tool wear diagnostics. The input to the system .are the 
multi-sensor raw data,and the cutting 4ool wear r measurem_ents collected during th^ experimental phase as 
well as the machining category. The term "machining category" represents the combinations of ( r^achming 

45 parameters that are .used 'for, collecting the training data.. When the machjning parameters vary considerably, 
the sensor features, also* vary widely., Hence, ..withpui. any. knowledge about..the. machining .conditions, and 
with just the sensory; data it- is difficult tc- establish .a diagnostic, mapping befweeathe sensory features £nd 
the cutting tool state. Subsequently, the knowledge, engineer divides the broad spectrum of machining" 
conditions into. severaj.% categories such that*, each, category represents, the, machining . parameters" that \ 

so resemble close to each other. An example, of machining -categories is. given in Tablet. . ; 
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Category 


Work piece 
material 


Ori 11 
size/type 


Speed 
1 " Crpm) 


Feed rate ^ 
(inch/rev.) 




Ductile* cast 
iron 


0.25" 
Tungsten 


1100-14Q0 


0 .006-0". 01 J 


r 


Ductile cast ' : 
" i ron 1 


6:25" - 

Tungsten 


1600-2000 


0.003-0^006 


2 


A3 6 steel 


0,25 a : High 
Speed Steel 


360G-39OQ 


Q.002-:0.004 

_ih : - 


3 


A>6 J ;steel 


L Speed Steel 


2360 T 2700 


0.002- 1 0.004 : 

j ' • * . ■ - 


4 l 


tA36 -steel 


0.25" High x 
Speed Steel 


■2300-2700., 


™ r. 11 - v---^r: 
0.008-0.010- 


5 


A36° steel t 


p.125 - High 
Speed Steel 


370Q-:;4ooo 1 


' 0.006-0 i 007 


6 ; * 


^A36*teel : 


:0.125" High 
Speed Steel 


3700-400a 


0.009-O.01O -| 


7 t 


A36 steel . v 


.0.12SV High \ . 
,Speed Steel 


29ig-iVgp 


p.oo ; 9^o r :o'ip- : '"'' I 



Table 1. Categorization of Widely Varying Hachitfing Cop 



The knowledge' base generation is divided into several sub-tasks: feature selection, -building the 
30 relational structure, determining feature state descriptors, and subjective estimation of^the conditional 
probability distributions. \ \_ 

1 . Feature Selection 

35 In Fig. 1, the spindle motor current sensor 114 is an AC current sensor. It is a noninvasive induction 
sensor that is simply clipped around the input power line. The 'AC current sensor output is rectified and low 
pass filtered in order to acquire low frequency energy of the sensor signal. This energy is directly 
proportional to the cutting torque exerted by the tool on the work piece. As the tool wears, the cutting force 
experienced by the tool increases. This increase in cutting force causes a corresponding increase inspindle 

40 motor current. The RMS value : pf spindle motor current thus becomes a valuable feature* for observing the 
trend in tool wear during the initial^hase Of drilling. . ' - : 1 - 

In addition to the RMS value, the change'' in RMS value compared^ a : fresh tool (ARMS) is : a useful 
feature as it indicates the temporal trend of the cutting torque! A faster rate of change' indicates a rapifrwear 
and hence shorter tool life." ; / . -..'*'.' ' " 

45 Another valuable source" of information about' the drilling processes the spindle" mounted strain' gage 
122. It indicates the thrust" force exerted by the tool" 100 on the Work" piece : 1 02. As the tooM 00 wears, the 
thrust forces increas^due to a constant feed rate and an increase in energy required to penetrate the work 
piece 102. Thd actual value used is the mearvbf trie strain gage'forc^ during the drilling of one-hole. 

Furthermore, as in the case* of the spindle motor current sensor 114; "the raw value (the mean) as weir 

so as its incremental change from a fresh 7 tool (Amearr) are used to predict the tool life. ■ * 1 0 ■ 

In the preferred embodiment the features extracted from the measurements of the spindle motor current 
sensor 11 4, and the strain gage sensor 122 are used to predict tool life and to continuously diagnose tool 
wear at the end of the predicted life of the tool. However, in an alternative embodiment other sensors,* for 
example the feed motor current sensor 110, the accelerometer 128, and th dynamometer 124, are used 

55 either in lieu of or in conjunction with either th spindle motor current sensor 114 or the strain gage s nsor 
.122. 

2. Building the relational structure 
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An Influence .. Diagram .Expert System (IDES) is an expert system in which . relationships between 
variables are quajjtativety and quantitatively represented and manipulated to arrive at spme" decision. An 
influence diagram, Lei, the relational structure, which is, a t graphical network 'and", theoretical structure for ; 
knowledge, provides probabilistic inference and expected value decisions, and has two hierarchical levels: , 
5 symbolic and quantitative. At the symbolic level, an influence diagram consists of . an acyclic, directed 
network with nodes representing relevant system state variables and , the arcs representing state variable" 

- interrelationships. For example., the diagram .in Fig. 5 is a simple influence diagram which schematically 
shows the probabilistic, influence, three variables, RMS, ARMS, and category, have on the state of a cutting 
tool (wearl). At the : quantitative level the mapping between "the relevant states are represented by discrete 

10 conditional probability density functions. For example, in the sample .influence^. diagram in Fig. 5, node 3 
would be quantitatively represented by the conditional probability { wearl \ RMS t ARMS, category }. Also, 
in a general model, the root nodes 1, 2, and 3 are quantitatively represented by their respective marginal 
probabilities. However, when the influence diagram technique is used in a real-time paradigm, they are 
represented by uniform distributions. A .diagnostic , inference, for example, the marginal probability of the 

75 cutting tool being worn, results from a.js^t of transformations applied to the diagram "in response to a 
specific query. The topological structure remaining after* the set of transformations representee answer to 
the query in symbolic form. At the quantitative level, these transformations would provide the answer to the 
query in numerical Jorm,. Three types of transformations used in this reduction process are Sensor Node 
Removal (table look-up operation), State Node Removal (applying summation rule of probability) and Arc 

20 Reversal (applying Bayes* rule). Once diagnostic inference has evaluated the likelihood of critical events, 
the best course, of action-must be determined. i<t: 
Fig. 6 is an influence, diagram for predicting tool life in the preferred embodiment. Verification uses an 
influence diagram with the same topology. However, at the numerical level, the state space of each feature 
and the subjective estimates of the probability distributions are different for different phases.. Hence.Jhe 

25 portion of the knowledge base corresponding to the influence diagram for the verification phase is different 
from the corresponding portion of the prediction phase. . 

The features "Arms and rms from the, spindle. motor current sensor influence the life 1 variable, an 
intermediate node in the influence diagram, "strain gage JEprqe features Amea/ii^and*yi7^an influence the Iife2 
variable, and ..machining parameters speed, feed, wprk-pi'ece^ material and* tool size influence the class 

30 variable. The latter is not a probabilistic influence but rather ^a^categprizati on. ~ An. example^ categorization 
based on these particular parameters is shown in Table 1 . How.ever,.other categorizations are possible. The 
class variable also influences the life! and Iife2 variables. The range'of values Jpr'/ifef 'ari6Jife2 is "long 
life" or "short life." Thus, given particular values~of each feature,' for example fms and Arms\ both values 
together indicate with a particular conditional probability .that the value of life 1 is -"long life." ; . ' /' z \ . 

35 Finally 1 life 1 and Hfe2 influence the //Ve variable, which, is an. output node in the influence diagram, the 
values available for the life Variable is "jong.life^and "short " life ".VSqth, of which are values' dependent on 
the type of tool in question. For example, "for 0.25" drills, based on experimental results, "Jong life" is 
defined as 100 holes. and "short life" as 10 holes^and.similariyjor Oil 25 h r drills, "long life" corresponds to 
50 holes and "short life "i to 7. hojes. However, in alternative embodiments these definitions may be different. 

40 , - ■ r ■ ' . * ; \ - r« i- . . 

- 3. Determining feature state descriptors 

The quantitative level of the influence diagram, as part of knowledge base 152, each node is divided 
into many states depending on the range of values, symbolic or numerical, the variable can take, The 

45 values that divide the range into many states .are referred to as.the state descriptors. The number of states 
varies. For example, for the life node 423,th v e lifel node 41*7/ and the Iife2 node 419, there are two states 
"long life" and "short life." The class node 421, has -as many states as the total number of categories, e.g. 
eight in the case of Table 1. In one embodiment, using the categories in Table 1, node 415 has "0.25-inch" 
and "0.125-inch" states, and node 413 has the states "steel" and "cast iron." Thus,, tlie, states for nodes 

so 423, 417, 419, 415 and 413, are^aU qualitative?The remaining nodes in^Fig: 6 have state descriptors defined 
by quantitative intervals. The states, for. nodes,. ^09 and 411 are obtained by dividing the honTial. operating 
ranges of speeds and feeds respectively, 1 - Wer "as" many intervals as needed to uniqueiy define each 
category in Table 1. For node 411, and Table 1, a state 0 represents feed rate less than 0.008 ipr, and a 
state 1 represents feed rates greater than or equal to 0.008 ipr. Similarly, in the case of node 409, a state 0 

55 is defined by drill speeds less than 3200 rpm, and a state 1 is defined as speeds greater than or equal to 
3200 rpm. In the case of feature nodes 401 through 407 the state descriptors are decided by analyzing and 
correlating the raw experimental data with the actual life for the tool for all the categories. 

Fig 7. is an illustration -of- raw. data collected for, twp,.category i> cases. It illustrates how to obtain the 

8 
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decision boundaries associated with the RMS feature. Thus, the figure represents a specific feature, RMS, 
of node 403, arid a specific category, category 5. The experimental data may be plotted in similar graphs 
for other features and categories: In Fig. 7, the RMS feature of spindle current signal for both drills is plotted 
as a function of number of holes drilled. Each point in the^graph represents the RMS value for signals 

5 collected forgone hole. The vertical dotted lines (labeled "Drill #1 Worn*' and "Drill #2 Worn", respectively) 
indicate the hole' number at which the drills' were actually declared as worn by wear measurement criterion. 
This hole will be referred to as the "wear hole." The circleslto the left of vertical line for drilf#i belong to 
the good tool condition and ones to the right belong to the worn tool condition for drill #1* Similarly 'drill #2 
is represented by plus symbols. Note that the life of Drill #1 is a "long life" whereas the life of Drill#2 is a 

70 "short life." The objective is to select one threshold (horizontal line called decision boundary) value that 
separates the "long life" state from the "short life" state for both the drills. For the prediction phase the 
average RMS value for the first three holes is used. Thus, the decision boundary for the prediction phase 
(labeled a) is selected so that the average RMS value for the first three holes will indicate whether a "lon f g 
life" is indicated or not. In the example of Fig. 7 the average RMS for the first three holes of Drill #1 is 

75 approximately l 0.45 and the average RMS for the first three "holes of Drill #2 is' approximately 0.36, therefore 
the selected decision boundary for RMS for Category #5 drills is 0.4. Similar "analysis is done for the 
verification phase. However, in the case of verification the, maximum feature values for the "drills up to the 
mid-life is used rather than the average for the first three holes. In the case of Drill #1 the max' value is 6.59, 
as indicated by the boxed ring. The max value for Drill #2 is 0.5. Therefore, the selected decision boundary 

20 for the verification phase for RMS for category #5 drills is 0.57. 

The decision boundaries for all features and all categories is done in the same manner as explained 
above for RMS and category #5. Thus, a quantitative knowledge base of decision boundaries is built. 



25 



45 



55 



4. Subjective estimation of the conditional probability distributions. 



The next step is to assign marginal probabilities associated with the decision boundaries. The marginal 
probability for .a given feature and category is a subjective evaluation on the part of the knowledge engineer. 
It indicates the probability that the decision boundary accurately predicts "long life." 

Having assigned the marginal" probabilities for all feature nodes 401 through 407, the knowledge" 
30 engineer determines conditional probabilities for the intermediate nodes 417 and 419. The class node, node 
421, representing machining category; is -assigned a binaryjdeterrhinistic conditional distribution (0 or 1) 
given any element of its coriditional cross product space consisting* of all the combinations of states of the 
predecessor nodes 409 through 41 {5 V c *' ' 

Using node 417 as an example, for 'each of the predecessoV features rms (node 403) and Arms (node 
35 401) a subjective weight oh a ? scale"of 0 to 1 is assigned based on how well each decision boundary 
separated the^ short life", from the "long life" states of dll the categories.' Let these weights' :j be w A 03 and 
iv 4 oi for rms and Arms, respectively. For eacH' category, critical' states We identified for "nodes'" 401 and 403, 
below which-me-feature-*-^^ which- the feature' values 

correspond to the "short life" state. For some category c\'let%*03 and s+01, be the critical states for nodes 
40 403 and 401 respectively. Let O403 and IHoi be the state spaces for nodes 403 and 401 respectively. The 
conditional probability assessment is made as follows. - 

For category c, the conditional probability, P, that the node 417 is in "short life" for (/ x J) « Q403 x fUoi 



P= — +A' 2 (Q)r * 'if i2s 403 ,j<s 401 y - - 

W 4Q3 W 401 ' ' 

,P= [ Wd01 +d 3 (QK' ■ l Xf 'i<s 40r Jzs 401} . ,../. 

s P=0.0; , if.i<s 403 ,j<s 40 t _ *\ 

where * . •_*-.. u. 

0 = (/ + J + S403 - s*o» +'1) ' '** '* ' ■'' : ' ' ' ••>*.■. : * 
and Ai, A2, and A3 are smoothing* Actors that 'gradually change the coriditional proBabilities in discrete' 
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steps as the features' states gradually change from 'those corresponding to a drill with a short life to, those, 
corresponding to -a drill with a long jife. They typically assame ; values ranging from 0;05\lo 0.1, and are* 
chosen subjectively*' All these subjective quantities are different for different categories :and this procedure; 
is repeated for all the Categories. While tuning the knowledge base, they rare iteratively: updated until -the.. 

5 system responds correctly for the training data. Typically no more than five iterations -:are. necessary, -The 
equations above are heuristic, and a person skilled in the art will realize other possible equations. 

As to the training phase with respect to tool wear diagnosis, Fig. 8 shows the topology of the influence 
diagram used by IDES 154 during the continuous diagnosis of drill wear in Step 617; of Fig. 3. A comparison 
with Fig. 6 indicates the^-similarity -between the influence diagram for tool life prediction and tool'.wear. 

70 diagnosis that takes' place at the- end oof the predicted life of the. tool. However,- interior, and -.output nodes 
differ. In the tool wear diagnosis influence diagram, the features derived from the spindle motor current 
influence the wearl variable and those derived from the strain gage force influence. therwear2 varjable. 
Those intermediate nodes influence the output node, wear. The values for wearl, wear2 and wear are 
"Fresh drill" and "Worn drill." Furthermore, at the numerical level, the decision boundaries and probabilities 

75 differ from the influence diagrams associated with tool life prediction. 

In the case of tool life prediction discussed above, the knowledge base is manually built by a knowledge 
engineer. The knowledge base for the continuous diagnosis of tool wear (of Step 621 of Fig. 3) is 
automatically created using a dynamic distance classifier. Fig. 9 shows the relationship of the dynamic 
- distance classifier to the total operation of the computer aided, manufacturing system haying an influence- 

20 diagram expert system. 

The. dynamic distance classifier generates the decision boundaries between the fresh and worn tool 
classes as well as the probability of correct classification for each feature. It uses .the Euclidian distance 
between the sample points and the mean of each class as the measure for classification. It bears the name 
"dynamic" due to an iterative update of the? fresh tool class'mean based; on "a new -separability -criterion 

25 defined below. The method , is generic ,in that r it . i s .applicable to drilling as well _as oiher machining 
operations, for example milling or turning. Furthermore, "applications outside of the computer aided 
manufacturing environment will be apparent to a person skilled in the art. 

The method assumes that the data is collected at regular intervals for a short period of time and the 
sample points, classified based on some tool wear criteria, represent a time history of the state of the tool. 

30 This assumption is not a restriction but is in-line with the true situation in a machining environment. 

1. Nomenclature 



The variables used in the method are defined as follows: 



35 


Co • 


Fresh cutting too! class 




Cl - 


Worn cutting tool class 




d- 


Decision boundary 




w - 


Actual sample # that divides classes 




T - 


Tolerance zone, in terms of number of samples, about w (e.g., T = {4,5,6,7} given w = 5) 


40 


Wmin " 


Minimum of ho r all ie T (e.g., iv mjn = 4 given T = {4,5,6,7}) 




X - 


Feature under evaluation for its classification capability 






Value of x at /" sample 




9i ' 


The feature mcai of class Q 




P- 


Probability of successful classification for the feature under evaluation 



45 

Separability Condition for Cutting Tool Classes 



Conventionally, in pattern classification using linear discriminant function F, two classes Co and Ci are 
separable if 
50 Ffa) < 0 for all x t € Co and 
F(x t ) > 0 for all x t e Ci , 

where X\ denotes the feature value for sample / and the decision boundary d is defined by 
F(d) = O. 

Since it is undesirable to have any fresh tool states to be classified as worn, the condition 
55 F(Xi) < 0 for all x x e Co 

is a requirement for correct classification. After w samples the tool is defined as worn based on some wear 
criterion. For conventional separability, the classifier diagnoses the tool to be worn for . each sample 
collected beyond the w* h sample. However, in case of robust diagnosis, at the first instance of wear 
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diagnosis, a decision would have been made; either to.- replace -the cutting tool or continue operation in 
degraded moder depending on" the economics- of that particular machining operation. Hence it is sufficient 
that there exists just one correct classification of the worn drill within a tolerance zone (T) ot a- few samples 
about the actual- worn tool sample w. The tolerance zone T may consist of a few. fresh tool samples before, 
the worn tool sample w (as indicated in the. definition of 7). Consequently, fresh and worn toolxlasses are. 
defined as separable if 'C^-t c ■ / ■ . •: 

1 - F(X\) < O for all Xj e Co and . r- ' >; . 1 

2. there exists U T such that F(xi}> 0,.i.e., x, € Ci •: - 

This relaxed but adequate separability. condition '.simplifies the search for a decision boundary _ that 
would provide robust-classification. For each individual feature, this .decision boundary. ' is, nothing^ but- a 
point. This one dimensional search for the decision boundary is considerably faster than the .search for the 
n-dimensional decision "surface where the features are;arranged into an nrdimensional vector. 



75 



3. Classifier Method -C o.--r - : - ?• 

The following pseudocode:; determines the decision boundary -and the marginal probability for, each 
feature. • -'^ -3. .- -■■ 



20 



For (x = *Feature( 1 ) ■ to Feature (ft) ; 



25 



Begin ?or loop 1 •"■ '. r 1 r 

( 1,) . Compute meam of each class (g,-; j =, 0,j;. :; 
yvr *'(2) : Jf (g 0 k '^^y^Sn^f^iufn (p ^ &) • antf exitV^ J 5 ; 
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(3) For (i = First sajnpl,e/ to Last sample;, . 
Do'classif ication using : decision rule, 
; - x, € C k *if \x if -g k \ = : min \x r .*-. : g-\', j = 0,1 
(A) Jf (x\ € C 0 for airi € T) then return (p = 0,;- and exit, 
(5) -If (x^ e C Q formal! i ■< w min ); <±nd 

(there exists i € T such that x, e C l ) then 



return d = ?°* 9x 



and exit. 



(6) While ( there exist i e T such that x x €"'C i y 

and (.there exists ; i ^ T such that x, C l ) 
BEGIN While " 



d 



9o + 9 X 
2 



P = 1 - — — # 

25 V min ~ 



^0 = 



. where s. = I S\ , 
:*• ; v S = {-^ ! 6 Ci for all i < w^} 

.* ? •*< . 

.. where. x t e C x for alJ i -< 



Jjf r?o - ?i/ then exit While loop. 
• " "Repeat step(~3) Only. . * c 

40 ; end ^ii ? / r ;; 'w^V ' V-.^- , ' 

(7) Pf= fx, ^ C 0 - for -321 i •< v -in ; and '( there exists i e T 
V_ -such th'at/'x* # C" .then' 

45 '•.-'-•••V retfuf^id- = 9o 4 ^ andp^ VT * 'and" exit J . • - 



. , . €t.i s e f e t u r h X4 ^ p j . ^ h'd " e x i t • ' 
END For loop 1 



Thus, the dynamic distance classifier computes decision boundaries (d) and their associated probabil- 
55 ities for correct clarification : (p) fon.eaqh:^ classifier relies 

on a predetermined - tolerance., zone;. "T, -for itirnel^vdLasnosis^ of tool. wear, and "the fact that samples ar 
sequential in time, to minimize the generation of false alarms. A false alarm is a diagnosis of a tool as.-being 
worn that is made while the tool -is still good,; .i.e., -before;, the itolerance zone. In step 6, by iteratively 
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updating the relevant parameters until the condition terminating the While loop has been satisfied, false 
alarms are avoided- Furthermore, p is weighted by s, the number of false alarms. * -7 : 

Returning to^Fig. 9, the heuristic rules, for computing the. conditional probabilities whe,n using the 
dynamic distance classifier differ from those described above with respect to tool life prediction and 
verification of too! life prediction. The' decision boundaries generated by the" dynamic classifier are assigned 
as the state boundaries for .each feature node. Using .the marginal probability of correct classification for 
each feature, some heuristic rules are used to automatically generate the'overall knowledge base for the 
real-time expert system. In the influence diagram of ; Fig. 8, for the rms node^803,-the decision boundaries 
for each category 0 to c, be denoted by a* ...*a n , where n £ c+1, and. the corresponding marginal 
probabilities for correct classification is denoted by fV r n (h can be less than c + 7, because the decision 
boundaries of some the categories may be made to coincide if they lie close to each other. Similarly, ... 
b n and n ... r n denote the decision boundaries and marginal probabilities for the Arms node 801. The node 
rms 803 is divided into n+ 1 states and the node Arms is. divided into m + 1 states. For example, node 803 
is in state 0 if x < ai , in state 1 if ai ^ x ^ a 2( ... in state n if x £ a n , where x is the rms feature value. The 
node 821 is in state 0 if the machining parameters correspond to category 0, ... . state c if the category is 
c. Also, the wearl node 817 is divided into two states, state* 0 for a fresh-tool and state 1 for a worn tqol. If 
k is the state of the category node, 821, the critical state corresponding to the decision boundary and 
probability for the rms feature are k 1 and r k1 , respectively, and fdr Ar/ns.'fc? and s^. 

Based on the above definitions, for each point (x x y) in the conditional cross-product space /Qsos x 
Osoi^) of nodes 803 and 801, where x and y denote respectively the current state of nodes 803 and 801, the 
following heuristic rules are used to compute the conditional probability that the wearl node is in state 7, 
i.e., worn tool state. 

Conditional probability that node wearl 817 is in "worn drill" state for (x x y) e (Q 8 p3 x Qboi) 



2S 



30 



35 



= 1.0 



= r 



k2 



Jo + JCj - x - y - 1 



m + n + 2 



172 



x -y 



m + n + 2 



if x i l^, y £ Jcj 

if x * k^'i y .<* k 2 

'if x y * k 2 

if x < k lt y < k 2 



40 



45 



so 



If the computed quantities using the above equations lie beyond [0,1], then trte corresponding limit 
values are substituted for the conditional probabilities. Also,, the* above formulation caa -be ^extended to 
include more than three nodes as the predecessors of the result nodes like wearl 817. ^ „ ... 

Returning to Fig. 9, the dynamic distance classifier computes the decision bound an esf and probabilities 
for correct classification. The Jarte< *are*-usedrin "the heuristic rales de.sjprjbed: "&bqye to generate^ the 
conditional probabilities. These, combined with the decision boundaries, and the relational structure, e.g. the 
influence diagram shown in Fig.8, are used to automatically 'generate the 1 fcnowfedge base. This step 
consists of writing the required information into the format for the knowledge base structure the Influence 
Diagram Expert System 154 expeqts^The, generated^nowledgelJbase j$, used according, to the method 
described above with respect to Figures 2 and 3. 

While thislnvention has been described with reference to illustrative embodiments, this description is 
not intended to be construed in a limiting sense. Various modifications and combinations of the illustrative 
embodiments, as well as other embodiments of the invention.^wil^ be. apparent to persons skilled in the art 
upon reference to the description. It is therefore intended that* the appended claims encorhfiasis any such 
modifications or embodiments. I «" - - ' 

In particular and in a first aspect of the present invention, a system for predicting the life of a cutting tool 
in an automated metal cutting machine, comprising: 



55 a plurality of sensors for Measuring* physical 'data* relating to a cuffing 'opferati'orv of -said machine;- and 

a computer conn cted'to said sensors,- corn prising a^featuf" 1 " xtracfor for transfoFmihoJ^aid physical 
data into feature values; ' l - ^ ' \L :V: * -"-"^ s <~ -:- ^Ict *o*'r y^-^- o; £>.:*.<"..■ ; . 
a tool life- predictor for arriving at a-'prediction of -the life-or : said cutting : tool i 'having : -"• - - .-* 
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an influence diagram, having an ..input level of feature nodes, intermediate^ nodes, connected to. said, 
feature nodes and to a machining class node, an output node connected to said intermediate nodes; 
a first module for. computing averages of feature values for each sensor; 

a second module responsive to said feature averages \ox making a classification of said feature average 
5 as indicative of a long life or a short life; and .,.,,,.'...'**" 

a third module connected to said influence diagram for combining .said classifications, to arrive at a 
classification and. a conditional probability at each intermediate, node and a classification at said output. node 
which is the predicted life of said cutting tool. „ _ .... 

In a further aspect of the present invention, a system includes a tool wear diagnostic module responsive 
10 to said predicted life of said cutting tool and to said.^verification of said predicted life pf.said cutting tool, for 
arriving at a diagnosis of wear of said cutting too!, having 

an influence diagram having an input level of feature nodes, intermediate nodes connected to said 
feature nodes and to a machining class node, an output node connected to said intermediate, nodes; 

a seventh module responsive to said feature values for making a classification of said feature values as' 
is indicative of a tool being worn or not worn; and 

a eighth module*connected to said influence diagram for combining said, classifications to, arrive at a 
classification and a conditional probability at each intermediate node and a classification at said" ojjtput node 
which is the diagnosis of wear of said cutting tool. 

20 Claims 

1. A system for predicting the life of a cutting tool in an automated metal cutting machine, comprising: 
a plurality of sensors for measuring physical data relating to a cutting operation, of said^ machine; and 

a computer connected to said sensors, comprising a feature extractor for transforming said physical 

25 data into feature values; - . -: r Ci ^ Vl • - , - * - r s p~ : \^ - nvo-: - v - A - ^ r- v - 

a tool life predictor responsive to said feature values to arrive at a predicted topi life of said .cutting "too I. 

2. The system of claim 1 and further comprising a tool life prediction verifier responsive to said predicted 
tool life and said feature values for arriving at a verification of saic predicted tool life. 

30 

3. The system of claim 1 or claim 2 and further comprising a tool wear diagnostic. module responsive to 
said predicted tool life, said verification and said feature .values to continuously diagnp.se said tool for a 
worn out condition to arrive at a diagnosis of tool wear of said cutting tool. 

35 4. The system of any preceding claim wherein said tool life, predictor, arrives at a prediction of the. life of 
said cutting tool, comprises: 

a knowledge base : haying.an influence diagram with an input level of feature nodes, intermediate nodes 
connected to said feature nodes and to a machining class node, and an output node connected to said 
intermediate nodes; ^ ... , . 

40 a first module for computing averages of feature values for each sensor; 

a second,, module, responsive to said- feature , averages -for making ri a. classification of ..said ..feature 
averages as indicative of a long life or a short life; and 

a third module connected to said influence. .diagram for combining said classifications^ to arrive at a, 
classification : and a conditional probability at each intermediate node and a classification at said' output 
45 node which is said predicted tool life of said cutting tool. 

5. The system of any of claims 2, 3 or 4. where said tool life, prediction verifier comprises: , . 
acknowledge base -having an influence, diagram, with" an input level of feature nodes,, intermediate" nodes 
connected to said feature nodes and to a machining class node, and an output node connected to 'said 

so intermediate nodes, ri.-c - * . . d : : ,- .- - . M ■ ; - 4 , , -„.. , „. 

a first module for computing maximum values of. said feature values for each sensor; , , . " f [ 

a second module responsive to said maximum values .for making a classification. oJL £9jd maximum 
value as indicative of a long Jife or. a. short -life; and-.- P . • ^ - - * - - 

a third module^ connected to ^said influence diagram fipj., combining said classlf^catiprjs to.^rrive at.a 

55 classification .and a conditional probabili^at .eaph intermediate. node. : and ( a. classification at said output 

node which is said- verification . 

6. The system of any of claims 3, 4 or 5 wherein said a tool wear diagnostic module comprises: 

14 
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a knowledge" base having an* influence diagram with an input level of feature nodes, intermediate-nodes 
connected to said feature nodes and to a machining class noderand an output node connected to said 
intermediate nodes; ' * ' ' * 

a first module ^responsive to said feature "values for making a classification of said feature value as 
s indicative of said tool being worn or not worn; and : 

a second module connected to said influence diagram for combining said classifications to arrive at a'' 
classification and a conditional probability at each intermediate node and a classification at said output 
node which is said diagnosis of tool wear for said cutting tool. ' * ; 

70 7. The system of any preceding claim, wherein said sensors comprise: ** t: 
a strain gage. * ... 

8. The system of any preceding claim, wherein said sensors comprise:' - \ - 
a spindle current sensor. • - - * 1 

75 ....... ... J 

9. The system of any of "preceding claim, wherein said sensors comprise: ■ ,L - 
a feed current sensor. * 

10. The system of any preceding claim, wherein said sensors comprise: 
20 a dynamometer. 

11. The system of any preceding claim, wherein said sensors comprise: 

an accelerometer for sensing vibrations in said machine: °' * ; , ' '• 

25 12. A method of operating a metal cutting machine having a cutting tool, comprising the steps -6f:- ! ^ * " 

(a) attaching at least one sensor to said metal cutting machine; " . • r " *' :%r : ~ *"~ 

(b) commencing a metal cutting procedure; 

(c) observing wearof said cutting tool and noting when said cutting tool is worn out according to 
some criterion; 

30 (d) collecting data from said sensors; 

(e) selecting from said sensors features which influence the wear of said cutting tool; and 

(f) using said selected features for predicting the iife b of said tool. ' ' 

13. The method of operating a metal cutting machine of Claim 12, further comprising the step of: 
35 (g) verifying the predicted life at the mid-point of the predicted life. 

14. The method of operating a metal cutting machine of Claim 12 or 13, wherein step (e) further comprises 
the step of: 1 .-.>■■ 

(e.1) determining decision boundaries for the data collected in step (d). - n 

15. The method of operating a metal cutting^ machine of Claim 14, ? where1ri step (e) further 'comprises' the 

steps of: --. ' *. f 

(e.2) selecting machining parameters for categorizing the rriaGhining operation; and * - '* : - 

(e.3) building an influence diagram for modelling the wear of said cutting tool' based on 'the data from 
45 said sensors. ' 

16. A method for ^learning in a real-time expert system, comprising the steps of: - 

(a) creating a relational structure corresponding to the influences that * first and a second variable 
have on a third variabfe; "' * * T ' ^ ' ' - ' ' * ' " 

so (b) obtaining data for said first variable.^ said second variable, and said third variable over some 

duration of time; ,^ - ' * * * . X . . l~ " ' ' 

(c) selecting a tolerance zon'e^around a sample that divides- classes; arid - : " * r: '- - ; ' ■ ■ 

(d) using a dynamic distance classifier, which is responsive" to 'said tdlerarice zone, to determine 
decision boundaries for ^"c? first and 'second variables Hand probabilities;- associated with said 

55 decision " boundaries, said' protiabilities ^iridicating' tH a s ' decision ^boundary 

classifies values of said first and second -variables with respect to' thervalue : of said third variable? 
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